First order phase transition with a logarithmic singularity 
in a model with absorbing states 
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Recently, Lipowski [cond-mat/0002378| investigated a stochastic lattice model which exhibits a 
discontinuous transition from an active phase into infinitely many absorbing states. Since the 
transition is accompanied by an apparent power-law singularity, it was conjectured that the model 
may combine features of first- and second-order phase transitions. In the present work it is shown 
that this singularity emerges as an artifact of the definition of the model in terms of products. 
Instead of a power law, we find a logarithmic singularity at the transition. Moreover, we generalize 
the model in such a way that the second-order phase transition becomes accessible. As expected, 
this transition belongs to the universality class of directed percolation. 



PACS numbers: 05.70.Ln, 64.60. Ak, 64.60.Ht 
I. INTRODUCTION 



In nonequilibrium statistical physics the study of phase 
transition continues to attract considerable attention [0 . 
In this context, continuous phase transition into absorb- 
ing states have been of particular interest ||^. It is be- 
lieved that absorbing-state transitions can be categorized 
into universality classes, the most prominent ones being 
directed percolation (DP) j|] and the so-called parity- 
conserving (PC) universality class [Q . On the other hand, 
various nonequilibrium models are known to exhibit a 
discontinuous phase transition . Especially in one spa- 
tial dimension first-order transitions require a very ro- 
bust mechanism in order to stabilize the ordered phases. 
As suggested in Ref. |^ , first-order transitions in one di- 
mension should be impossible under certain generic con- 
ditions if one of the ordered phases fiuctuates. 

Recently, Lipowski introduced a model with infinitely 
many absorbing states which exhibits a first-order tran- 
sition from a fluctuating active state into an absorbing 
phase 0J^. Remarkably, this transition takes place even 
in one spatial dimension. Unlike previously investigated 
models, the dynamic rules involve products of real- valued 
local variables. Surprisingly, the first-order transition is 
accompanied by an apparent power-law singularity of the 
stationary particle density, suggesting that the model 
may combine features of continuous and discontinuous 
phase transitions. This observation collides with the 
commonly accepted belief that there are no power-law 
singularities in first-order phase transitions. The aim of 
this work is to study the origin of this unusual type of 
singularity in more detail. 

The model considered in Ref. is defined as follows. 
Each site i of a given lattice is connected to n neighboring 
sites j €< i >. Each bond carries a real- valued variable 
Wij = Wji G (—1/2, +1/2). The set of all bond variables 
specifies the state of the system. A site is considered to 
be active if the product of all adjacent bond variables is 
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FIG. 1. Lattice geometries used in Ref. [8]. Left: 
One-dimensional triangular ladder, interpreted as a linear 
chain with next-to- nearest neighbor interactions. Right: 
Two-dimensional square lattice. 



smaller than a certain control parameter r: 



n 



< r . 



(1) 



Initially all bond variables are uniformly distributed be- 
tween — 1/2 and +1/2. The model evolves by random 
sequential updates according to the following dynamic 
rules. For each update attempt a site is randomly se- 
lected. If it is active, all adjacent bond variables are 
replaced by new random numbers distributed between 
— 1/2 and +1/2. Lipowski considered the case of n = 4 
neighbors using a one-dimensional triangular ladder and 
a two-dimensional square lattice (see Fig. |l|). 

Due to the use of real- valued local variables, the model 
has infinitely many absorbing states. Moreover, the dy- 
namic rules are invariant under certain gauge transforma- 
tions. For example, we may invert the sign of all bond 
variables along a closed contour without changing the 
pattern of activity in a given configuration. The impli- 
cations of this type of gauge invariance are not yet fully 
understood. 

In order to understand the existence of a phase transi- 
tion, let us consider two extremal situations. On the one 
hand, for r > 2^" it is obvious that all sites are active 
during the entire temporal evolution. On the other hand, 
for r < there is a finite probability to generate bond 
variables with \wij \ < 2"~^|r|. This means that certain 
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pairs of sites i and j remain inactive forever, irrespective 
of the values of the other bond variables. Thus, the pro- 
cess continuously " switches off' certain pairs of sites and 
therefore approaches an absorbing configuration within 
an exponentially short time. 

Interestingly, for r=0 the model is still in the active 
phase with a non-vanishing stationary density of active 
sites po > (see Fig. ||). Thus, the spreading process un- 
dergoes a discontinuous phase transition at r = 0. Even 
more surprisingly, the stationary density Ps{r) does not 
decrease linearly versus as r — > 0, instead the slope of 
the curve seems to diverge. This observation led Lipowski 
to the conjecture that the model may combine features 
of discontinuous and continuous transitions, calling for a 
power-law behavior of the form 

Ps{r) ^ pQ + ar^ for r > , 

Ps{r) =0 for r < 0. ^ ' 

where a is a certain factor and /3 is the critical ex- 
ponent associated with the order parameter p. Per- 
forming numerical simulations Lipowski found the esti- 
mates pq ~ 0.314827, 13 ~ 0.66(3) in one dimension and 
Pq = 0.358, /3 = 0.58(1) in two dimensions, respectively. 
Moreover, he observed that the dynamical critical expo- 
nent z — Vi_/ V]i^ ~ 0.2 is very small. 

In the present paper we propose a different explanation 
of the diverging slope in Fig. ^ based on the following 
arguments: 

1. The singularity of the slope in Fig. |^ emerges as 
an artifact of the definition of activity in terms 
of a product. While the bond variables are uni- 
formly distributed, the probability distribution of 
the product diverges for r — > 0, leading to a singu- 
larity of the slope. An explicit formula is derived, 
showing that the slope diverges logarithmically as 
r ^ 0. In particular, there there is no power law 
of the form (|). 

2. Redefining the control parameter, the model dis- 
plays a conventional first-order phase transition 
without a diverging slope. Moreover, there is no 
diverging length scale at the transition. 

3. Since for r < the system is immediately driven 
towards one of the absorbing states, the continu- 
ous transition may be thought of as being hidden 
in the inaccessible region r < 0. In order to support 
this point of view, we generalize the model in such 
a way that the continuous transition is shifted to 
the accessible region r > 0. As expected, the tran- 
sition belongs to the universality class of directed 
percolation. 
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FIG. 2. Stationary density ps as a function of the param- 
eter r in one (ID) and two (2D) spatial dimensions. The inset 
zooms the region where the curves terminate. The terminal 
points at p(0) = po are marked by the bold dots. 

II. LOGARITHMIC SINGULARITY 

In this Section we demonstrate that the singularity of 
the slope in Fig. || is a consequence of the multiplica- 
tive definition of activity in Eq. (|l|). To this end we 
consider the reactivation probability W{r) that a site re- 
mains active after an update. In a spreading process, 
the reactivation probability provides a good measure of 
the effective spreading rate by which nearest neighbors 
will be activated. In most spreading processes near the 
transition, the reactivation probability varies to lowest 
order linearly with the control parameter of the model. 
Therefore, a power law of the form ps ~ {W{r) — W{rc))^ 
immediately implies the same power law in terms of the 
control parameter ps {r — rc)^ ■ In the present model, 
however, W{r) is a nonlinear function at the transition. 
Therefore, it does matter whether r or W{r) is used as 
control parameter. 

Our explanation relies on the assumption that W{r) is 
the "true" control parameter of the model. In terms of 
W{r) the model exhibits a regular first order phase tran- 
sition without singularity, i.e., the density varies linearly 
with W(r) — W{rc)- The apparent singularity in terms 
of r originates solely in the non-analytic behavior of the 
function W{r) in the limit r 0. 

In order to compute W{r) let us consider n indepen- 
dent random numbers zi,Z2,... , Zn drawn from a flat 
distribution between —1/2 and -f 1/2. Clearly, the prob- 
ability P^^^{z)dz to find one of these random numbers 
between z and z + dz is constant for \z\ < 1/2. How- 
ever, the product z = Y[i=i of the random numbers is 
not uniformly distributed. The probability distribution 
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p("^ (z) can be computed iteratively by 

p('=+i)(z)= dz' dz" P'^''\z')5{z- z'z") 

, 1 

' dz' ^P'^''\z/z')Q{z' -\z\2^) 



= 2 



rf2'ip('=)(|z|/z') 



(3) 



with P^^\z) = 1, leading to the exact result 
(n — 1)! 

Thus, the probability W{r) to reactivate an updated site 
is given by 



W{r) 



dzP^^\z) 



-+r( 8-8(log, 16|r| 



(5) 



4(log,16|r|f--(log, 16| 



For small values of r, we therefore expect the stationary 
density to be given by 



p,(r) ~ po + a{w{t)-]^ 
Ps{r) = 



for r > , 
for r < . 



(6) 



where A is a fit parameter. This formula replaces the 
power law in Eq. (^), implying that the slope of Psif) 
diverges logarithmically as — (logg 16|r|)^ for r ^ in 
any dimension. This explains why Lipowski's results in 
one and two dimensions are so similar. 

In order to support the validity of Eq. (|), we per- 
formed Monte Carlo simulations (see Table D). Our es- 
timate for the stationary density po in one dimension 
deviates slightly from the value quoted in This de- 
viation may be explained as follows. On the one hand, 
the random number generator plays a crucial role. Since 
most algorithms generate internally an integer random 
number, the output is often quantized in steps of about 
10~®, leading to wrong results if r is very small. On the 
other hand, the machine precision itself limits the accu- 
racy. This problem can be avoided by storing log(|it;ij|) 
and sgn(u'y ) instead of Wij and turning the product in 



Eq. (|l|) into a sum of logarithms. Taking these technical 
subtleties into account, we obtain a different estimate. 

As shown in Fig. ^, the results for the stationary den- 
sity in one and two dimensions are very similar. More- 
over, for small values of r the curvature of the lines is in 
fair agreement with the theoretical prediction of Eq. 
In any case, the possibility of a power-law singularity can 
be ruled out. As expected, there is no singularity if ps is 
plotted against W{r) - W{Q) (see Fig. |). 





d= 1 


d = 2 


number of sites 


10"^ 


1000^ 


simulation time 


10* . . . 10'= 


10* . . . 10" 


stationary density po 


0.31512(3) 


0.35905(2) 


fit parameter A 


1.4(1) 


1.1(1) 



TABLE I. Numerical estimates. 
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FIG. 3. Log-log plot of the stationary density Ps(r) as 
a function of the parameter r, compared to the (vertically 
shifted) function W{r) ~ 1/2. The dotted straight line visu- 
alizes the failure of the power- law conjecture proposed in [8]. 
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FIG. 4. Stationary density ps as a function of W{r) — 1/2. 
As can be seen, there is no singularity near the transition. 
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FIG. 5. Phase diagram of the generalized model. The 
solid (dashed) line indicates the phase transition line in one 
(two) dimensions. 



IV. CONCLUSIONS 

The common feature of the models introduced by 
Lipowski et al. is the use of products of real- valued local 
variables in the definition of the dynamic rules. The mod- 
els display very interesting phenomena which are usually 
not observed in models with linear local rules. As we have 
shown in the present paper, the use of nonlinear rules is 
crucial since nonlinear functions of random numbers may 
not be uniformly distributed. A special situation emerges 
if the distribution exhibits a singularity. In this case a 
tiny change of the control parameter may lead to a dra- 
matic variation of the order parameter. For the model 
investigated in jsj , such a singularity is responsible for the 
diverging slope in Fig. Thus, in contrast to a previous 
conjecture, the model does not combine features of first- 
and second-order phase transitions. In particular, the 
model is not critical and the correlation length remains 
finite as r — > 0. We expect that similar phenomena may 
emerge whenever the dynamic rules are defined in terms 
of nonlinear functions of real- valued random variables. 



III. CONTINUOUS TRANSITION IN A 
GENERALIZED MODEL 



As outlined in the Introduction, the first-order phase 
transition at r = is induced by frozen pairs of inactive 
sites. Since for any negative value of r there is a finite 
probability to generate such pairs during the temporal 
evolution, the system is driven into one of the absorb- 
ing states within an exponentially short time. Roughly 
speaking, the spreading process is switched off as soon as 
r < 0. Thus, the continuous transition, which is expected 
to exist in ordinary spreading processes, may be thought 
of as being hidden in the inaccessible region r < 0. 

In order to support this point of view, we generalize 
the model in such a way that the continuous transition is 
shifted to the accessible region r > 0. This can be done 
by introducing two control parameters ri , r2 and consid- 
ering a site to be active if the product of adjacent bond 
variables lies in the interval {—ri, r2)- Clearly, this model 
includes the original one as a special case. Moreover, the 
phase diagram is symmetric under exchange ri ^ r2. In 
order to avoid frozen pairs of sites, we will assume that 
both parameters are positive. Obviously, for very small 
values of ri and r2 the spreading probability is so small 
that the model will be in the absorbing phase. Increasing 
ri and r2, we observe a continuous transition from the 
absorbing to the active phase, as shown in Fig. ||. 

We verified that the critical behavior along the entire 
transition line belongs to the universality class of directed 
percolation. Note that the Z2 symmetry along the diag- 
onal ri = r2 does not lead a different type of transition 
since it is not a symmetry of the order parameter. 
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